The Rotating- Wave Approximation: Consistency and Applicability from an 

Quantum System Analysis 



Open 



Chris Fleming, 1 '* N. I. Cummings, 1 '^ Charis Anastopoulos, 2 ' * and B. L. Hu 1 ' 3 '^ 

1 Joint Quantum Institute and Department of Physics, 
University of Maryland, College Park, Maryland 20742-4111, USA 
Department of Physics, University of Patras, 26500 Patras, Greece 
' Maryland Center for Fundamental Physics, University of Maryland, College Park, Maryland 20742-4111, USA 

(Dated: March 8, 2010) 

We provide an in-depth and thorough treatment of the validity of the rotating-wave approximation 
(RWA) in an open quantum system. We find that when it is introduced after tracing out the 
environment, all timescales of the open system are correctly reproduced, but the details of the 
quantum state may not be. The RWA made before the trace is more problematic: it results in 
incorrect values for environmentally-induced shifts to system frequencies, and the resulting theory 
I has no Markovian limit. We point out that great care must be taken when coupling two open 

f-^ systems together under the RWA. Though the RWA can yield a master equation of Lindblad form 

^vq similar to what one might get in the Markovian limit with white noise, the master equation for 

the two coupled systems is not a simple combination of the master equation for each system, as is 
possible in the Markovian limit. Such a naive combination yields inaccurate dynamics. To obtain 
the correct master equation for the composite system a proper consideration of the non-Markovian 
<^ dynamics is required. 

00 

Q_l I. INTRODUCTION 

+-» 

A. Background 

a 

^■i The rotating-wave approximation (RWA) is used in many places in the study of open quantum systems, particularly 
in the field of quantum optics (see for example [1—3]), but the validity of the approximation is treated in depth far less 
often. There are actually two distinct rotating-wave approximations both in widespread use : 1) the 'pre-trace' (preT) 
RWA, which consists of modifying the interaction Hamiltonian by dropping the so-called counter-rotating terms that 
are quickly oscillating in the Dirac picture; and 2) the 'post-trace' (posT) RWA, which is obtained by neglecting terms 
in the master equation for the reduced density matrix that are quickly oscillating in the Dirac picture (see, e.g., [2, 3] 
C and [1, 6] respectively). Agarwal has carried out a systematic study [4-6] differentiating between these two RWAs 

and addressed their validity for atom-field interactions and spontaneous emission processes. More recently, various 
(f-) authors have claimed some features of the RWA that may limit its applicability, which we will now discuss. 

The most widely acknowledged problem with the preT RWA seems to be that it yields incorrect frequency shifts 
in the atomic energy levels, so that it is not suitable for calculating environmentally induced level shifts or induced 
cooperative frequency shifts [4-6]. West and Lindenberg [8] found that the reduced system dynamics obtained from 
J> the pre-trace RWA does not have a Markovian limit 1 . Finally, Ford and O'Connell [7] have raised concerns that in 
general the total Hamiltonian obtained by the preT RWA does not have a spectrum which is bounded below, and 
they suggest that this limits the applicability of the approximation to first-order transition amplitudes. 

Other authors have raised a very different sort of concern about the preT-RWA Hamiltonian for coupling of a 
localized system to a quantized field, that it may produce spurious causality violation in the calculations. Consider, 
for example, a two-level atom in the dipole approximation interacting with the electromagnetic field. The multipolar 
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1 These authors also claim that there is no fluctuation-dissipation theorem for this model, a statement that the present authors cannot 
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form of the Hamiltonian for this interaction is 

H = d.D(R)=J2f:^ l (d-% s ) Ke-* + CT _e-*) 0^) (1.1) 

in the Dirac picture, and in the preT RWA the terms with frequency u) + would be neglected. However, with these 

terms dropped the interaction can no longer be expressed in terms of the local field variable D(R) [ ]. Indeed, a 
numerical study of a three-atom problem [14] found that noncausal terms appear when the preT RWA is used, unless 
one makes the ad hoc modification of extending frequency integrals to —00. The preT RWA may then misrepresent 
the effects of retarded propagation in the electromagnetic field, which suggests problems with causality in the study 
of multipartite systems. 

Moreover, the Glauber detector model [10], long used in photodetection theory and quantum optics, uses the preT 
RWA, and it might give rise to quantum correlations between spacelike separated events that do not represent the 
effects of actual entanglement. The effective status of preT RWA in Glauber's theory is debated: some authors 
have shown that photodetection probabilities at short times appear to violate causality [11], and modifications to 
Glauber's photodetection theory have been suggested [11, 12], while others indicate that a different form of the RWA 
in photodetection theory can guarantee causality [1 3]. Our interest in this problem partly arose from finding how the 
imposition of preT RWA affects the range of validity of results from the calculation of the entanglement dynamics of 
two atoms interacting with a common quantum field at large atomic separation [15]. 

Some form of the RWA is also invoked in discussions of open quantum systems in the context of quantum optics 
and atom-field interactions. Indeed, the standard derivation of the Born-Markov master equation in general requires 
an RWA to render it in Lindblad form [ 1 ] , thus providing a completely-positive dynamical map (for all states at all 
times) as is conveniently assumed in many quantum information theory discussions. 

With these two distinct RWAs in widespread use while some open questions remain about their limitations and 
fallacies, we find it useful to carry out a systematic analysis of the consistency and applicability of the RWAs in the 
modern language of open quantum systems. This would be fitting in view of the fact that researchers today have to 
tackle problems beyond those of level population and dissipation rates to deal with more subtle issues such as quantum 
decoherence and entanglement dynamics and perform more demanding tasks such as quantum state tomography and 
engineering. 

B. Main Results 

This major findings in this work are: 

1. The RWA may be sufficient or insufficient depending on what information is desired about the system. For the 
perturbative relaxation rates either the pre-trace or post-trace RWA is sufficient. To obtain the environmentally 
induced shifts in system frequencies, only the post-trace RWA is sufficient. In order to get more detailed 
information about the evolution of the quantum state and the asymptotic steady state neither RWA is sufficient 
in general. 

2. Even in a system where the underlying environmental noise is colored, the RWA can yield a Lindblad-form 
master equation similar to what one might get in the Markovian limit with white noise. However, when two 
such systems are coupled together, one cannot obtain the correct dynamics by naively adding together the 
dissipative terms of the master equations, as would be possible in the Markovian limit. While the master 
equation obtained by such a naive addition has a mathematically valid form, it does not generally yield the 
dynamics of the system obtained from a microscopic derivation of the master equation for the combined system 
(even after the RWA). There are sufficiently simple systems and bath correlations for which a less judicious 
application of the RWA can nevertheless produce an adequate master equation. However, this is no guarantee 
for validity in more complex situations. 

3. The pre-trace RWA does not in general have a Markovian limit. 

This paper is organized as follows: In Sec. II we review the RWA as it applies to interactions within a closed system. 
In Sec. Ill we discuss the consistency and applicability of the post-trace RWA. In Sec. IV we discuss the same for the 
pre-trace RWA. In the last section we conclude with some discussions. 
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II. THE ROTATING- WAVE APPROXIMATION IN CLOSED SYSTEMS 

We first examine the rotating-wave approximation as ordinarily applied to a closed system consisting of several 
interacting subsystems. Its wave function evolves as 



dt 



m = 



'■HI 



(11.1) 



under the total Hamiltonian H = Ho+Hi where Ho represents the sum of all uncoupled subsystems and Hi represents 
the part from the subsystems coupling. One seeks to solve the eigenvalue problem 



H \oj) = uj \ui) , 

for eigenstates \uj) with eigenvalues w by perturbing off the free eigensystem 
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where ui^ are the eigenvalues of the uncoupled system. The non- degenerate first-order corrections are then 
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and for the degenerate corrections one must find the correct linear combination of degenerate states |/), which exist 



solely in the degenerate subspace 
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This kind of analysis should also extend to nearly- degenerate subspaces where the basis corrections in Eq. (II. 6) would 
nearly diverge. 

In the rotating-wave approximation one only considers components of the interaction Hi which oscillate least 
rapidly in the interaction picture. If these terms are stationary, e.g. at resonance, then they will include the correct 
non-degenerate first-order frequencies, Eq. (II. 5), as well as the correct characteristic equation, Eq. (II. 8), which 
determines the degenerate first-order frequencies and zeroth-order energy states. What the RWA generally neglects 
are the first-order basis corrections, Eq. (II. 6). If the RWA terms are non-stationary, then they will include the most 
nearly-degenerate first-order frequencies and also neglect their first order basis corrections. Therefore the RWA has 
limited correspondence to perturbation theory as long as all terms that are close to resonance have been retained 
in the interaction Hamiltonian. If, however, the RWA is made such that there are neglected terms that are near 
resonance, then the correspondence fails even in the weak coupling regime (because the missing first-order basis 
corrections become large, and even the first-order eigenvalue corrections are inaccurate). Moreover, for subsystems 
with a multiplicity of timescales, there may be several near-resonance frequencies for which a proper application of 
RWA would have to take into account. 

In the rest of this paper we consider open quantum systems and divide our attention between two cases: the post- 
and pre- trace RWA. 



III. THE POST-TRACE ROTATING-WAVE APPROXIMATION 



An open-system is produced by coarse graining over the environmental degrees of freedom and its dynamics is 
described by a master equation governing the reduced density matrix (RDM) p. In the case that the effect of the 



4 



environment can be represented by a quantum white noise 2 the master equation has time-local coefficients and 
assumes a Lindblad form: 

p = -i[K + V,p}+V{p}, (III.1) 
V{p} = V nm U n pel - X - {ele n , P }) , (III.2) 

where H is the Hamiltonian of the free system and V is a correction introduced by the environment, not necessarily 
to be renormalized in its entirety as it may contain nontrivial features such as diffusion components. ~D is the 
"dissipator" super-operator: T> nm is a positive-definite and Hermitian coefficient matrix and e n denotes a particular 
basis of representation for the dissipator. 

More generally, one can motivate the Lindblad equation [21, 22] from Choi's theorem on completely-positive maps 
j 10]. The Lindblad equation classifies all completely-positive master equations without time dependence in their 
coefficients. Following the generalization of Choi's theorem on merely Hermitian-preserving maps, one can quickly 
prove that any Hermitian and trace-preserving master equation of analogous time-local form must have a pseudo- 
Lindblad form with merely Hermitian T>. 



A. General Features 



The master equation for the RDM of an open quantum system can be written in the form 

p = C {p} + 5C{p} , (III.3) 

where 

£ {p} = -*[H,p] (III.4) 

is the term that would already be present in the unitary dynamics of the closed system and SC represents the 
dissipative corrections introduced by coupling to the environment. Such open-system master equations can be derived 
perturbatively, in the system-environment interaction, with a variety of different techniques [23-25] . The rotating- wave 
approximation is often introduced, not in a purely perturbatively derived master equation, but in one derived via the 
Born-Markov approximation. To second order in the system-environment interaction, the Born-Markov approximation 
is consistent with weak coupling perturbation, even well outside of the Markovian limit. Therefore, given that the 
RWA will only be applicable to second order, it is of no consequence if one starts from the Born-Markov approximation 
or a more rigorous perturbative analysis. One only has to keep in mind that RWA has no reliance upon any kind of 
Markovian approximation. 

The post-trace (posT) RWA effectively consists of only considering the parts of the super-operator SC which 
commute with the free system propagation super-operator 

G (t){p} = e- lHt pe +,Ht . (IIL5) 

If we consider evaluating our master equation coefficients in the energy basis, which in pseudo-Lindblad form amounts 
to resolving 

T^ki-ji = Vnm e„ \<Ji) e m (III.6) 

nm 

with \ujj) representing the system energy eigenstate with frequency Uj, then the RWA essentially amounts to projecting 
out the diagonal of this Hermitian matrix, i.e. terms with ujk — ^i — or equivalently ioi —loj = — u;/. For the 

perturbative master equation to second order in the system-environment coupling, which is assumed to be weak, these 
diagonal entries will settle to positive values and therefore this projection yields a master equation of the Lindblad 
form. We will refer to the master equation obtained this way as the RWA-Lindblad equation. 

Such a Lindblad projection is only reasonable because the system-environment coupling is assumed to be weak and 
the projection is performed in the energy basis. As we will presently show, the eigen-operators of the Liouvillian have 



2 Some authors refer to an environment or its noise as Markovian or non-Markovian, but we maintain that the terminology 'Markovian' 
or 'non-Markovian' be saved for the description of stochastic processes, not the nature of noise. 
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the form \wj) (u>k\ plus corrections at second order in the coupling strength so that the discrepancy introduced by 
the dropped terms is small for sufficiently weak coupling. The RWA-Lindblad equation is not fully equivalent to the 
weak coupling master equation, but it generates an evolution which is very close to that of the weak coupling master 
equation in a perturbative sense. 



1. Correspondence with Perturbation Theory 

Assuming the system-environment coupling to be weak, one can solve an open-system master equation in a pertur- 
bative fashion. One seeks a solution to the eigenvalue problem 

Ca = ftr, (III.7) 

for the matrices cr which would evolve with characteristic rate /, given the perturbative expansions 

C = Cq + SC, (111.8) 

<Tij = \uji){uj.j | + 8cr ij H , (III. 9) 

fij = -ioj ij + 6f ij + -- - , (III. 10) 

where the zeroth order terms are correctly set to match the free evolution of the system and Wjj = uji~ujj. For Gaussian 
noise, which has very well behaved asymptotic properties, the perturbative corrections are at minimum second order 
in the system-environment interaction. Therefore we will speak of all lowest order perturbative corrections as being 
second order. The second-order constraint upon our eigenvalue problem is then 

(ujk\ 8C{\ui)(ujj\} \uji) = - Lou) (uk\ S&ij \oji) + SfijSij-ki ■ (III.11) 

Evaluating the components of this equation yields the non-degenerate corrections. 

Ser^ \oji) = - — — J " ' where uj 1} ^ uj h , (111.12) 

—l\0Jij — Ul k i) 

s/ij = (wiiarflwiXwiDK-) • (HI.13) 

For era and fa the system is degenerate, and one must solve the characteristic equation (Pauli master equation) 

W<x = Sfa, (III. 14) 

(wi|W|w 3 -> = M^K-XwiDIwi) , (111-15) 
<Tj = <r \uji) , (111.16) 

for the branching under the perturbation. 

The essential point is that the perturbative corrections to the eigenvalues are entirely captured by the post-trace 
RWA, while the perturbative basis corrections are entirely neglected. The lack of basis perturbation can be readily seen 
in the thermal state of the system. The system evolving under the RWA-Lindblad equation will relax into the thermal 
state described by the Boltzmann density matrix p cx H . But this is generally not how most systems would actually 
thcrmalize in a Hamiltonian formulation: Given the full system + bath and interaction Hamiltonian, H + Hb + Hi, 
the system is expected to relax into a state described by the reduced density matrix p cx (e~P( H+HB+Hl ^ B , assuming 
the bath to be very large 3 . This state typically only reduces to Boltzmann form in the limit of zero system-bath 
interaction strength; however, an infinitesimal interaction strength would imply an infinite relaxation time, so we 
cannot speak self-consistently of the system relaxing to the thermal state the RWA-Lindblad equation predicts. 

This property suggests a limitation to the applicability of preT RWA in the study of entanglement dynamics. In 
a multipartite system with components interacting only through the bath, this precludes the presence of asymptotic 
residual entanglement. This is in contrast to the asymptotic behavior of bipartite systems in quantum Brownian 
motion-see, for example [29, 30]. 



That the entire system + environment appears thermalized in its reduction is proven to second order in Ref. [28]. 
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2. RWA Fails When Perturbation Theory Fails 



The second-order master equation should be valid when the second-order corrections provided by 8C are small as 
compared to the unperturbed dynamics generated by Cq. Let us denote the strength of the dissipative corrections 
generically by the frequency 7d, then the weak coupling condition is 



as uiij corresponds to the eigenvalues of Co (not Ui which corresponds to the eigenvalues of the unperturbed Schrodingcr 
equation) . 

The RWA-Lindblad equation does not directly correspond to the second-order master equation, but more correctly 
to the second-order solutions of said master equation. For the second-order solutions to be valid the coupling must 
not only be weak in the above sense but also in the following sense 



which justifies the perturbative solutions. One cannot have near degeneracy in the energy level splittings or the naive 
perturbative solutions, which the RWA-Lindblad equation corresponds to, will fail. Perfect degeneracy is acceptable; 
the RWA-Lindblad equation retains these terms in the Pauli master equation for instance. But near-degeneracy needs 
to be treated in a manner analogous to degeneracy; the nearly degenerate subspace should be diagonalized. With 
the second-order master equation this is still possible, but with the RWA-Lindblad equation these terms have been 
discarded and one is left with an invalid master equation. Thus, there can be situations where the weak-coupling 
condition is satisfied while the posT RWA condition is not. 

This problem arises, for example, in cavity QED. A two-level atom of frequency f2 coupled to a resonant intracavity 
field mode will result in an energy spectrum of the composite system that has the form of the harmonic oscillator 
with each level split in two by the Rabi frequency ^JnQyR of the dressed states, 



If the intracavity field is coupled to the field outside the cavity, this becomes an open quantum system. If the 
intracavity field is coupled weakly enough to the atom, then the system will be in the weak-coupling regime of cavity 
QED and the vacuum Rabi frequency ilvii will be small compared to 7d ■ In this case the posT RWA procedure does 
not, strictly speaking, apply as was noted by Scala et. al. [16]. 

In such a case one can still do a partial RWA, neglecting terms that oscillate much faster than 7d and keeping those 
that are slower. This still leaves the master equation in pseudo-Lindblad form; however, assuming these timescales 
are sufficiently slow and the spectrum of environmental noise is sufficiently flat, one may be able to make an effective 
Markovian approximation for the remaining pseudo-Lindblad terms (even if one might not have been valid for original 
master equation due to the faster system dynamics that have been ignored in the posT RWA) to recover a Lindblad- 
form master equation. Scala et. al. argued this is the case for cavity QED with a low-temperature bath [17]. 



Before discussing the non-Markovian properties of the RWA-Lindblad equation, it is beneficial to get a clearer 
definition of non-Markovian, as the term is applied in a variety of incompatible ways in the mathematics, physics and 
chemistry literature. In an open quantum system perspective, the environment acts as a stochastic source or noise, 
of both quantum and classical (thermal) origins, influencing the system. This noise, which imparts fluctuations and 
dissipation, can be multiplicative when there is nonlinear interaction amongst the constituents in the environment 
and colored when there is temporal correlation (memory) which reflects the dynamical timescales of the environment 
(see discussions e.g., in [26, 33]). When the effects of a thermal reservoir with quantum mechanical degrees of freedom 
may be represented by a quantum white noise source, the open quantum system's noise kernel and damping kernel 
(the dissipation kernel's anti-derivative) will have delta function correlations. Noise which induces local dissipation 
will only necessarily have delta correlations for the damping kernel. Non-Markovian refers to physical processes with 
memories. In an open-system framework the influence of the coarse-grained environment can engender non-local 
noise correlations even with local dissipation or in the high temperature regime. White noise only enters at high 
temperature and with local dissipation. 

In the Markovian limit the timescales of the environment are taken to be much shorter than the timescales of 
the system. Thus in an open-system perspective one cannot simplistically refer to an environment as Markovian or 




(111.17) 



7d < min \cjij - ujm 



(111.18) 




3. Non-Markovian processes, nonlocal dissipation, colored noise 
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non-Markovian in isolation from the dynamics of the system which it influences but is measured in reference to. A 
noise may appear 'Markovian' only because its characteristic time scales cannot be resolved by the system. Some 
authors use the terminology 'Markovian noise' or 'non-Markovian noise' to describe the nature of the noise. We 
prefer to refer the nature or properties of noise simply as white or colored while reserve the terminology Markovian 
or non-Markovian for the description of processes, without or with memory, as it should be. 

Another common habit is the use of Markovian in reference to the master equation itself, if it is time-local and espe- 
cially if it is additionally time-homogeneous. This should more correctly be referred to as a Markovian representation 
or simply avoided. Markovian processes produce Markovian representations, but not all Markovian representations 
arise from Markovian processes. A master equation which arises from a non-Markovian process, even if in Markovian 
representation, is not sufficient to generate the dynamics of multi-time correlations. By definition, Markovian pro- 
cesses must give rise to a Quantum Regression Theorem in which higher order correlations can be determined from 
lower order correlations and ultimately everything can be resolved via the master equation. But with a non-Markovian 
process there are corrections to the Quantum Regression Theorem[18, 19]. 



A master equation for the reduced density matrix is time-local or convolutionless if the time derivative of the 
(reduced) density matrix at a time t is expressed only in terms of the (reduced) density matrix at that time p(t) (as 
opposed to an intcgro-diffcrential equation depending on the past history of p) . A master equation can be time-local 
and yet non-Markovian. A well-known example is the HPZ master equation [26] where the coefficients are obtained 
from solutions of ordinary differential equations, the non-Markovian features manifest in the nonlocal dissipation 
arising from the back-action of the environment with colored noise on the system. 

The corresponding quantum Langevin equations (Heisenberg equations of motion for system operators with the 
reservoir operators integrated out) will remain non-local in time. The class of time-local master equations is distin- 
guished from the class of master equations with constant coefficients such as the familiar Lindblad master equation 
on the one hand, and the time-nonlocal (integro-differential) equations on the other hand. The latter form is often 
encountered in a projection operator formalism with little or no coarse-graining of the environment. (See, e.g., [35]). 

A non-Markovian master equation will exhibit memory in a manner that is not naively apparent in its time-local 
representation. This is true even of the RWA-Lindblad equation, despite the fact that it is of Lindblad form as one 
would have in a Markovian limit. Let us consider the master equations for the reduced density matrices of the (open) 
systems A and B with Hamiltonians Ha,Hb each coupled to a dissipative environment of equivalent influence. We 
then have the open system 



where the Hamiltonians are those of the free systems and the corrections to the Liouville operator are introduced via 
interaction with the dissipative environment. In the Hamiltonian formalism one can simply add two Hamiltonians and 
arrive at another Hamiltonian, though one might be motivated to fix the energy spectrum through renormalization. 
One cannot do this with non-Markovian Liouville operators. E.g. given some coupling Hab between subsystems one 
cannot simply add dissipative terms. 



The above (incorrect) master equation is in general completely different from the correct open-system master equa- 
tion derived from first principles. For non-Markovian processes, the environmental contributions have a nontrivial 
dependence (due to memory effects) upon the system's dynamics through their couplings. If one changes the system 
Hamiltonian then one must also change the environmental contributions to be compatible with the history these new 
terms will create. This is how memory exhibits itself in a time-local representation. Moreover, one must also take 
into account whether or not the dissipative environments are separate or shared. If the dissipative environment is 
shared then the two subsystems can interact via environmental back-reaction. This effect is also missed when simply 
combining the Liouville operators. 

In a general non-Markovian master equation, the problems of the above incorrect master equation would be readily 
apparent as it would likely violate positivity, uncertainty, etc. Positivity violation will not occur when adding RWA- 
Lindblad terms, but the mistake has only become more subtle and therefore more dangerous. The master equation 
might be completely positive, but it does not correspond to the dynamics of the physical system considered. This is 
the key point. 



4- Constant-in-time, Time-local and Time-nonlocal Master Equations 



Pa = -«[Ha,Pa] + S£ -a{pa}> 
Pb = -»[H b ,Pb] + 5£ B {p B }, 



(III. 19) 
(111.20) 



Pa+b ^ -z[H A +H B + H A b, Pa+b]+<*£a{Pa+b} + <5£b{Pa+b}> 



(111.21) 
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This issue has already been commented on in the context of cavity QED. The often-used master equation includes 
the Hamiltonian for the atom, intracavity field, and atom-field interaction, but the dissipator used is exactly that 
of an empty cavity with dissipation plus that of an atom spontaneously emitting into empty space, so that the 
situation is just that depicted in Eq. (III. 21). And, indeed, if one begins instead with the atom-cavity system and 
derives the microscopic master equation using the standard technique [1], one finds that the master equation has a 
different dissipative term [16, 17]. As explained in [17] if the spectrum of environmental noise is sufficiently flat then 
the difference is suppressed, which explains the success of the standard cavity QED master equation. But not so 
otherwise, which is something often overlooked. 



B. Application to the Two-Level Atom 



Here we consider a two- level system with <j z Hamiltonian and energy level splitting f2, bilinear ly coupled to a 
thermal reservoir via a cr x coupling. This would, for example, model a two-level atom coupled to the electromagnetic 
field in the dipole approximation. Denoting our reduced density matrix 



P-+ P— 

one can compute the second-order master equation [28] and place it into the form 

r 



d_ 
dt 

d_ 
dt 



P++ 
P— 

P+- 

P-+ 



cosh( J^) 



— e 1 2t 



-e 2t 

n_ 

-e 2t 



P++ 
P— 



-r - i(n - sn) 



p+- 
p-+ 



with decoherence rate (here also the half thermalization rate) and energy level shift 

r = r(0), 



SQ 
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7T 



poo 




/ deV 




la 



(111.22) 

(111.23) 
(111.24) 

(111.25) 
(111.26) 



in terms of the phenomenological decoherence rate function r(w). V denotes the Cauchy principal value which 
regulates contained poles from contributing to the integral. 

This master equation and those that follow are exact to second order, only the coefficients have been allowed to 
relax to their asymptotic values. The relaxation occurs quickly, within the system and bath timescales, as compared 
to their effect, which occurs in the coupling timescalc. Therefore, when considering properly correlated initial states 
which do not jolt, it is safe to consider this "late-time" regime. 

In terms of the microscopically derived damping kernel 7(w), the anti-derivative of the dissipation kernel, the 
decoherence rate can be expressed 



r(fi) = 7 (ft)ftcoth 



2T 



(111.27) 



Regardless of system-environment coupling, the damping kernel is effectively constant for Ohmic coupling, which 
along with high temperature is responsible for thermal white noise. For linear coupling to the collective positions of a 
bath of harmonic oscillators the dissipation kernel has no more temperature dependence than the system-environment 
coupling itself. 

The post-trace RWA here amounts to neglecting the dynamical interaction between p+- and p-+. To second order 
in the coupling, the only effect of this is to neglect a perturbative amount of phase information pertaining to their 
damped oscillations, i.e. the perturbative change of basis. The asymptotic state works out to be exactly the same in 
cither case. Thus under the specific conditions leading to these results the post-trace RWA can be viewed as largely 
acceptable and somewhat innocuous. 



C. Application to Quantum Brownian Motion 



The exact Quantum Brownian Motion (QBM) master equation [2G] for the reduced dynamics of an oscillator 
bilincarly coupled (position-position) to a thermal reservoir of harmonic oscillators can be written 

p = -*[H R , p] - tT[x, {p, p}} - M D pp [ X , [x, p}} - D xp [ X , [p, p]] , (111.28) 
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where Hr denotes the system Hamiltonian with renormalized frequency coefficient, T the dissipation coefficient, D pp 
and D xp the regular and anomalous diffusion coefficients. This master equation will have a stationary limit if the 
noise correlation is not excessively widespread in time, e.g. a regulated ohmic coupling is perfectly suitable. To lowest 
order in the coupling, the late-time expressions for these coefficients can be determined from the weak coupling master 
equation [28] to be 



r 

Or 



r(fi), 



= ci- 



D rm 



D, 



deV 



ft 2 



r ( ())0<-ot.h( 2 *M 



deV 



1 



^(s) e coth^ 



£ 

2T 



(111.29) 
(111.30) 

(111.31) 

(111.32) 



in terms of the phenomenological dissipation function, which is proportional to the dissipation kernel at second order. 
In QBM with linear coupling to the bath, the dissipation coefficient will always have no more temperature dependence 
than the system-environment coupling. Ohmic coupling here will imply a constant dissipation function, though a cutoff 
is required for this theory. The frequency shift will have a linear cutoff sensitivity which can be renormalized, but the 
anomalous diffusion D xp has a logarithmic cutoff sensitivity. 

For this problem, the Fokker-Planck equation for the probability distribution function (Wigner function) W in a 
phase space representation presents a much cleaner picture with simple solutions [ ]. 



W 

q 
v 



{v T Hq+ V T DV} W, 

[x,p) T , 
d_ _9 1 T 
dx ' dp] ' 



(111.33) 
(111.34) 

(111.35) 



The matrices H and D are the homogeneous and diffusion coefficient matrices respectively. 



n 

D 



o -if 

Mfi| 2r 



— -D 
' 2 M Dpp 



(111.36) 
(111.37) 



Do not confuse the homogeneous generator with the Hamiltonian; they differ by some frequency rcnormalization and 
the dissipation T. From hereon we will assume the system frequency to be properly renormalized such that S7r = fi 
in the stationary limit. The diffusion matrix contains two components: the regular diffusion D pp and an anomalous 
anti-diffusion D xp which keeps the position uncertainty insensitive to high frequency. 

If one expresses the QBM master equation in terms of ladder operators, the pseudo-Lindblad coefficient matrix can 
be calculated to be 



1 

O 



Dpp iD X p Dpp 



(111.38) 



The rotating-wave approximation then constitutes projecting out the diagonal of this matrix, which will be positive 
definite. Transforming back into the phase space representation, the Fokker-Plank coefficients become 



'Hrwa 
Drwa 



r — i- 



2MQ? 







MD P 
2^ 



(IIL39) 
(111.40) 



The anomalous diffusion coefficient vanishes entirely while the dissipation and regular diffusion coefficients are both 
broken in half, with the missing half reappearing as an analogous coefficient of the master equation. 

The role of the homogeneous coefficients are to generate the homogeneous propagator e~ m . The RWA homogeneous 
coefficients are just slightly off in both the oscillation rates and phase; the dissipation rates are entirely correct. 
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Compare the characteristic frequencies of the two matrices 

h = T±i V / n 2 -r 2 , (IH.41) 

h RWA = T±iQ. (111.42) 

The diffusion coefficients are relatively more mangled given that the anomalous coefficient is entirely absent; despite 
what is amiss in the early calculations of the weak coupling QBM master equation, this coefficient does exist at lowest 
order in the system-environment coupling. The effect of diffusion is only present in the second cumulant or covariance 
of the Wigner function. For this stationary master equation, the evolution of the covariance is simply 



ait) 



-fH 



cr(0) - cr(oo)] « 



-tu 1 



+ er(oo) 



where the stationary covariance is determined by the Lyapunov equation 

Uer(oo) +cr{oo)n T = 2D. 
We can easily compare the stationary covariances. 



er(co) 
crwa(oo) 



(±D -D ) 

Mil 2 \2T 1J PP ^xp) 





IT W 



i i n 

MSI 2 2Y 1J VV 





M 







(111.43) 
(111.44) 

(111.45) 
(111.46) 



Amazingly the only difference in the stationary state will come from the lack of an anomalous diffusion coefficient. 
This contribution will ultimately be lower order in the coupling, due to the T _1 prefactor before D pp , and therefore 
its absence is acceptable perturbatively. 



IV. THE PRE-TRACE ROTATING- WAVE APPROXIMATION 



A. General Features 



1. Inconsistency of Approximation 



Let us consider a bilinear interaction Hamiltonian Hi between a system observable L and the collective environment 
observable 1. 

Hj = LI, (IV.l) 

For each of these operators, assuming them to be completely non-stationary, there is a gross raising and lowering 
decomposition L± given by 



L + = J2^ L ^ 



such that 



i<3 



L = L++L _, 
J = i(L+-L_) ; 
L, 



>jjj 



i4 



(IV.2) 
(IV.3) 



(IV.4) 
(IV.5) 
(IV.6) 



where L and J will be two relevant observables. For position coupling with a harmonic oscillator the decomposition 
becomes 



L = x . 



L+ = 



1 



L = 



J = 



V2MQ 
1 

V2MQ 



a, 



MQ, 



P: 



(IV.7) 
(IV.8) 

(IV.9) 

(IV.10) 
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L = 


0~x ; 




1 


>± = 




J = 


-<Ty 



and for er x coupling with a <x 2 Hamiltonian (two-level system) we have 

(IV.11) 
(IV.12) 
(IV.13) 

Now consider coupling the system to an environment made of a large number of harmonic oscillators in their 
collective positions. Let us furthermore assume the system coupling is like that of the above harmonic oscillator or 
two-level system such that it is characterized by a single frequency uj. 

1 = 5> fe x fc , (IV. 14) 

k 

Hx = ^== { ( L + ^ + L 4) + (L+ 4 + L a fe ) } , (IV.15) 

where xj. is the environment position operator with ladder operator a. k , energy e k and mass m^. In the interaction 
picture we have the interaction Hamiltonian 

H T (i) = ]T ° fc (L+ a fc e +«(-^)* + L aj e -'C«-«*>*) + 



V ; Cfc f L+ a f fe e +*(*>+^)* + L _ a fe e -*("+**)^ . (IV.16) 

An often utilized pre-trace rotating-wave approximation (preT RWA) is to neglect the second terms, conventionally 
referred to as counter-rotating terms, as they are deemed more rapidly oscillating than the first. However, this is only 
true in a mode-by-mode comparison. Keeping terms of frequency \u> — 2uj\ = uj while discarding terms of frequency 
\oj + 0| = uj serves no good purpose. There is no a priori sense in which this is an approximation at all, unless the 
only environment modes which exist are near resonance. 

A true bandwidth approximation which does what the preT RWA claims would instead modify the interaction 
Hamiltonian of Eq. (IV.16) (before tracing out the environment) by neglecting all the "rapid" terms that oscillate 
with a frequency outside some frequency band Auj in the interaction picture while retaining all the slower terms . 
The resulting Hamiltonian would be 



Hl (f) = y / ° k (l+ a fc e -K"-s*)« + L_ at e -»("-«»)A + 



Y , ° k ( L+ at e+^+^ t + L a fc er^^) (IV.17) 



£fc=0 
Auj— uj 



Note that if Aw < u>, then the bandwidth approximated Hamiltonian would have no counter- rotating terms . Fur- 
thermore, if the environment were such that all environmental frequencies s k lie in a band around resonance with 
\e k — oj\ < uj, then a bandwidth approximation using this band would be equivalent to dropping all counter-rotating 
terms. However, in the general case the two approximations are inequivalent, and simply dropping all counter-rotating 
terms is inconsistent. 

It is also important to note that if the bandwidth approximation of Eq. (IV.17) is performed with Aw chosen such 
that all near-degenerate terms are retained, then this is just the sort of RWA we discussed in Sect. II. The only 
difference is that the environment is to be traced out at the end of the calculation. However, such a bandwidth 
approximation would render the problem more difficult to solve than simply calculating a full perturbative solution. 



4 Note that this "bandwidth" Hamiltonian does not arise from restricting the field to some bandwidth of modes around the resonance 
frequency. 
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2. Noise and the Markovian Limit 



The Hamiltonian obtained after RWA is not generally an approximation of the full interaction Hamiltonian for 
reservoirs. It is nonetheless a linear Hamiltonian interaction with a thermal reservoir and will affect dissipation, 
decoherence, thermalization, etc. Therefore it still possesses some of the same character to the original model. 

Back in terms of observables, the RWA interaction Hamiltonian takes the form 

H : = jLVcfcX. + JjVc^, (IV.18) 

2 2 ^ m k^k 

k k 

and thus it describes a different but related set of system variables coupled to a different, but related set of bath 
variables. This results in two quantum noise sources 

Irwa = ^^CfcX fc , (IV.19) 
Jrwa - ly,Ck- —, (IV.20) 

which have not only autocorrelations in and of themselves but cross-correlations between themselves. Perhaps more 
clearly, if we consider the original damping kernel with one L-coupled source 



1 




(IV.21) 



7(w) = 7(0;) 

then the RWA damping kernel with both L and J-coupled sources becomes 

(IV.22) 



/ \ 7M 
7rwaM = — 



1 —i sign (us ) 
-zsign(w) 1 



with reference to the original damping kernel "/(uj). The damping kernel, here in Fourier space, is defined (as in 
Ref. [28]) as the anti-derivative of the multivariate dissipation kernel, which is itself the imaginary part of the mul- 
tivariate noise correlation in the time domain. The diagonal components come from the self-correlations (l(t) 1(t))b 
and (j(£)j(y))B, while the off-diagonal components come from the cross-correlations (l(t)j(r))B and (j(i)l(r))B. 

There is a subtle pathology in the cross-correlations of these two noise sources. The RWA interaction is an example 
of couplings to different kinds of bath observables with strong cross-coupling. Such couplings do always not admit a 
Markovian limit. The reason for this is because in addition to high temperature, the white noise limit also requires 
a local damping kernel, i.e. one constant in the Fourier domain. This is not a problem with one noise source as one 
can typically choose an appropriate coupling, e.g. Ohmic, such that the damping kernel will work out to be local. 
But with multivariate noise one must make all components of the damping tensor local, including new kinds of terms 
which arise from the cross-correlations. Whether or not this is possible depends in part upon any relation between 
the self-correlations and the cross-correlations. 

For the RWA damping tensor, if we make the diagonal components local with what was Ohmic coupling, then the 
off-diagonal components will appear highly non-local like sign(w). But if we were to choose a coupling as to make the 
off-diagonal components local, then the diagonal components will necessarily be highly non-local. There is no choice 
of coupling which can give us white noise. This problem with the white noise limit of pre-trace RWA has been noted 
before [8]. 

3. Correspondence with Perturbation Theory 

The perturbative correspondence between the pre-trace RWA and the original model is a bit more complicated to 
demonstrate. Let us start with the second-order corrections for our simple separable coupling without any sort of 
RWA following Ref. [28]. 

{<jj k \5C{\u)i){ujj\}\u)i) = (uj k \ L [A(uijei) + A(u)ij)] (u!j \ L \wi) 

h 

-Shi L K> A M (u> h \L \ui) , (IV.23) 

h 
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with late-time master equation coefficients 



A(w) 
«(*) 



1 , « 
2 aM - 2^ 

(1(*)1(0)> B , 



deV 



(IV.24) 
(IV.25) 



where a is the quantum noise correlation for our stationary bath. These corrections capture all of the second-order 
relaxation rates, perturbative frequency shifts, and basis corrections. 

As we have discussed, the post-trace RWA essentially considers taking only the diagonal entries where Wki = u>ij, 
and under appropriate conditions this is sufficient to reproduce all of the perturbative frequency shifts and relaxation 
rates but not the basis corrections. 



(uji + u)\ SC{\LUi)(ujj\} \tjOj + u) = (uji + ui\ L \u)i) 2 Re[A(w)] (ujj + uj\ L \ujj) 



■E 



A(a/)|(^+cy|Lk>r 



(IV.26) 



One can see that the first terms, which directly correspond to the pseudo-Lindblad dissipator, are now only determined 
by the real part of A(w) or the characteristic function of the noise correlation a(ui). This function is always positive 
by Bochner's theorem. 

The pre-trace RWA master equation has four related sets of terms because of the two correlated noise sources. But 
as far as these diagonal terms are concerned, which determine the perturbative timescales, one can essentially consider 
a master equation of the same form but with the modified coefficients 



. . , 1_, . i f + °°, l + sign(u;) sign(e) , 
Arwa(w) = - a(w)-- J de , 



1 



a(e) . 



(IV.27) 



The real part, which determines the relaxation rates, remains unchanged. But the imaginary part, which determines 
the energy level shifts, is very different. So while the post-trace RWA can correctly produce all of the perturbative 
timescales, the pre-trace RWA can only produce the relaxation rates, consistent with what has been found in earlier 
specific cases [4, 6]. 



4- Non-Markovian Nature of the Master Equation 



It is necessary to point out that, although the pre-trace RWA can often produce a master equation of Lindblad 
form, the coefficients are inherently non-Markovian. Even though the master equation is in a convolutionless form, 
the coefficients themselves contain integrals over the system's history alongside nonlocal correlations of the noise. As 
such, they cannot be universally applied to different systems (which would have different histories) even if one only 
wants the relaxation rates. This was easy to notice for the post-trace RWA as the correctly derived master equation 
coefficients would come out to be completely different. Here the reason is much the same. 

For instance, let us consider an oscillator system with x coupling to the environment. The accuracy of the pre-trace 
RWA decay rates stems from a correct , a raising and lowering operator decomposition of x. This leads to a different 
but related model of environmental interaction with x and p-coupled noise. We have proven that the perturbative 
decay timescales will work out to be equivalent, but only by using the raising and lowering properties. If we couple 
this oscillator to additional degrees of freedom in some larger system, then and a are no longer ensured to be raising 
and lowering operators for the new energy eigenstates of the system. Once this criterion has been broken, the proof 
fails to apply and all coefficients of the misapplied master equation will likely be wrong. A correct pre-trace RWA 
interaction would have to involve a raising and lowering operator decomposition which utilizes the full Hamiltonian 
of the larger system. 
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B. Application to the Two-Level Atom 

Utilizing the second-order master equation, we find the RWA interaction Hamiltonian yields 



P+- 
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(IV.28) 
(IV.29) 

(IV.30) 



In addition to differing from the frequency shift without the RWA, it also contains a higher order cutoff sensitivity. 
For approximately local dissipation, the sensitivity was logarithmic but is now linear. 



C. Application to Quantum Brownian Motion 



Again utilizing the second-order master equation, we find the RWA interaction Hamiltonian yields 
P = -*[H^,p]-*r[x,{p > p}] -MD pp [x,[x,p\] -D* p [x,[p,p\], 
with new frequency shift and anomalous diffusion coefficient 



% = 0- 



D* = 

xp 



1 



deV 
deP 



e 

e-Q 



r( £ ). 



e-n 



r(e) e coth 



-) 
2TJ 



(IV.31) 

(IV.32) 
(IV.33) 



For what was Ohmic coupling in the original model, the frequency shift has a different but still linear cutoff sensitivity. 
However the anomalous diffusion coefficient now also has a linear cutoff sensitivity. If the cutoff is very large, this 
could be very problematic. 



D. A Multipartite Example 

Let us say that we have an array of, otherwise non-interacting, parallel q-bits all with cr z Hamiltonians. A simple 
dipole interaction can be represented with the bilinear interaction Hamiltonian 

Hi = ^<r Xn l„ (IV.34) 

n 

where cr Xn is the x spin component of the n th qubit and 1„ is its corresponding collective environment coupling. The 
environmental coupling for a qubit at location r„ is 

In = ]T g k {e +ik r " + e"^ a k } , (IV.35) 

k 

with gk oc for an electro-magnetic field environment such as discussed in Ref. [2, 34]. The resultant damping 

kernel corresponding to the (l n (t) l m (r))B correlation is 

7nm(w) = 70 sinc(r„ m o;) , (IV.36) 

where r nm = r n — r m and therefore the damping is Ohmic or local for the autocorrelations where r nn = 0. The 
cross-correlations, which are strictly nonlocal, vanish in the limit of large distance separation. 
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The pre-trace RWA interaction, which was considered for two qubits in Ref. [15], introduces the duplication of 
quantum noise sources 



RWA 



1 



jRWA 
•in 



k 



e +lk r " a k - e- ik r " a k 



(IV.37) 
(IV.38) 



and one must now consider the correlations between all such operators. The resultant damping kernels can be 
organized 



7nmM 
4 



1 

-i sign(w) 



-j sign(cj) 
1 



(IV.39) 



with reference to the original damping kernel j nm (uj). The scenario is much the same. The damping rates will be 
correct, while the frequency shifts and basis corrections (including asymptotic entanglement) will be incorrect. There 
is no longer a white noise limit, even when the qubits are distantly separated. 

However this remains a fairly reasonable physical theory, as the RWA interaction itself was fairly reasonable. The 
cross-correlations between different qubits still vanishes for large separations. The second-order master equation, 
being determined by the second-order cumulants or two-time correlations, will therefore reduce to that of qubits 
coupled to independent environments in the large separation limit. 



V. DISCUSSION 



In this work we have systematically examined the rotating-wave approximation by using a master equation for 
an open quantum system weakly- coupled to a general environment. There are, in fact, two distinct rotating-wave 
approximations: The pre-trace RWA is an approximation performed on the interaction Hamiltonian before the en- 
vironment is traced out which yields a somewhat modified Hamiltonian dynamics from which the reduced dynamics 
can be derived. The post-trace RWA is performed on the master equation for the reduced density matrix after the 
environment has already been traced out. In [28] we have derived the master equation for a system with a discrete 
energy spectrum weakly coupled to a stationary quantum bath. Using this general framework we have compared the 
master equation describing the dynamics of this open system under the pre and post-trace RWA to the dynamics 
without making such an approximation. We have specifically addressed the master equations for a two-level system 
and for a linear oscillator, two models in which the RWA is often invoked. 

We find the post-trace RWA to be more innocuous than the pre-trace RWA. It can be seen as an often reasonable 
approximation in which the full weak-coupling Liouvillian (which is time-local and of pseudo-Lindblad form) is pro- 
jected onto a Lindblad-form Liouvillian. We call the resulting master equation the RWA-Lindblad equation. We find 
that for a general open quantum system the post-trace RWA will yield exactly the same timescales as perturbative 
solutions of the weak-coupling master equation. The perturbative corrections to eigen-operators of the Liouvillian are 
neglected in the RWA-Lindblad equation, so the predicted quantum state will differ. In particular, the steady state 
solution of the RWA-Lindblad master equation will differ from the true steady state by an amount that is perturbative 
in the coupling. These results are consistent with what Agarwal found [4-6] for the two-level atom and the linear 
oscillator in the Born-Markov approximation. One context in which the discrepancy in the steady state could be 
important is examining the late-time behavior of entanglement dynamics at low temperature. When the system is 
bipartite and the ground state is separable, the RWA-Lindblad equation will give an asymptotically separable state, 
whereas the weak-coupling master equation leaves open the possibility of asymptotic entanglement. In view of our 
findings, we can say that generally the post-trace RWA is suitable if one only wants the perturbative timescales of 
the dynamics, but it may not be appropriate if one wants more detailed information about the quantum state of 
the system, as it misses some corrections introduced by the coupling to the environment, and it will also not be 
appropriate when perturbative timescales fail, i.e. for near-resonance in the energy level splittings. 

We find the pre-trace RWA to be more problematic. When the environment contains many frequencies with a 
spread comparable to the frequencies of the system, then the pre-trace RWA in general does not provide a faithful 
representation of the true solution. We also find that the pre-trace RWA results in two strongly correlated sources of 
environmental noise that together have no Markovian limit. The cross-correlations between the noise sources are such 
that if the autocorrelations are white then the cross-correlations are strongly colored. This issue has been noticed 
before [8]. Finally, we have shown that, unlike the post-trace RWA, the pre-trace RWA in general does not correctly 
obtain all perturbative timescales for the dynamics, yielding incorrect frequency shifts. This finding based on a more 
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extended theoretical framework agrees with results obtained for specific cases studied before for the two-level atom 
[ ]■ 

We caution that the way the RWA is applied also matters. For Markovian processes in closed and open systems, 
certain liberties can be taken with the master equation; terms can be mixed and matched rather haphazardly. A 
Markovian process will generally produce a master equation of Markovian representation and Lindblad form. The 
RWA can also produce a master equation of Markovian representation and Lindblad form, and therefore one might 
assume that these liberties can also be taken with the RWA-Lindblad master equation. But this is not the case as the 
underlying stochastic process remains non-Markovian and the master equation coefficients contain memory despite 
their Markovian representation. A haphazard construction of RWA-Lindblad master equations for multipartite systems 
can produce an evolution which is completely positive and yet totally unphysical. 

The pre-trace RWA as applied to quantum Brownian motion is additionally interesting because the interaction 
Hamiltonian in this case coincides with the Gardiner-Collett Hamiltonian [31] used to model the coupling between 
the intracavity field of a high-finesse electromagnetic cavity and external field modes. However, in this case the form 
of the Hamiltonian is not the result of a rotating-wave approximation and can be derived [32] from a "modes-of-thc- 
universe" approach for a cavity with a partially transmitting mirror in the limit that transmission is weak. So in this 
case the pre-trace RWA Hamiltonian corresponds to the physical Hamiltonian of an actual system, and the solutions 
of the master equation have relevance directly, rather than as an approximation. 

A single two-level atom is clearly a particularly simple quantum system. As such, some of the shortcomings of the 
RWA are obscured in this case. We find that the post-trace RWA gives the correct equilibrium state for a thermal 
environment in this case, in addition to the correct timescales, though it does miss some of the corrections to the 
transient quantum evolution that can be obtained from the weak-coupling master equation without the RWA. Thus, 
if one's theoretical investigations are limited to those features that it captures correctly, then the post-trace RWA 
may be a suitable approach. 

Our general analysis provides a reasonably good characterization of the validity of the RWA for constructing a master 
equation; however, many questions remain about its effects on the internal dynamics of more complex systems, such 
as the effect of either form of RWA on causality in multi-atom models with spatial separation [14] and its effect on 
entanglement dynamics [15]. 
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